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SECTION  I 


INTRODUCTION 


BACKGROUND 

The  objective  of  this  report  is  to  present  a  set  of  "general"  stress 
intensity  factor  solutions  for  several  cracked  hole  configurations.  As 
described  later,  the  crack  surfaces  are  loaded  with  an  arbitrary  crack 
face  pressure  in  a  manner  which  allows  the  results  to  be  applied  to  many 
other  practical  loading  conditions.  Solutions  are  given  for  the  through- 
the-thickness,  corner,  and  surface  cracked  hole  geometries  shown  schemat¬ 
ically  in  Figure  1. 

It  is  commonly  known  that  such  cracks  often  originate  at  fastener 
holes  in  mechanically  fastened  joints,  and  may  represent  a  failure  source 
for  aircraft  structures.  The  United  States  Air  Force,  for  example,  has 
adopted  fracture  mechanics  based  design  criteria  which  require  new  aircraft 
to  resist  failure  from  pre-existent  cracks  assumed  present  at  fastener 
holes  (1) .  Such  analyses  require  stress  intensity  factor  (K^)  solutions 
for  the  assumed  flaw  geometry.  Since  the  fastener  hole  can  be  subjected 
to  a  variety  of  complex  loadings,  many  different  stress  intensity  factor 
solutions  may  be  required.  In  addition  to  varying  degrees  of  load  transfer 
through  the  fastener,  for  example,  residual  stresses  may  be  introduced 
around  the  hole  (by  coldworking  or  by  use  of  interference  fit  fasteners) 
to  improve  the  joint  fatigue  life  (2-4).  The  loading  conditions  may  be 
further  complicated  by  bending,  thermal,  or  dynamic  effects. 

This  report  describes  a  series  of  "crack  face  pressure"  stress  intensity 
factor  solutions  which  may  be  readily  applied  to  a  wide  variety  of  such 
loading  conditions.  The  remainder  of  this  section  overviews  the  superposi¬ 
tion  techn  ,ue  which  :ovides  the  general  utility  for  the  crack  face  pressure 
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solutions  presented  here.  Subsequent  sections  describe  the  analysis  techniques 
used  to  obtain  the  crack  face  pressure  solutions,  discuss  the  results,  and 
demonstrate  their  application  through  several  example  problems. 


LINEAR  SUPERPOSITION 


In  order  to  describe  the  nature  of  the  crack  face  pressure  stress  intensity 
factors,  and  demonstrate  their  generality,  the  well  known  linear  superposi¬ 
tion  procedure  shown  schematically  in  Figures  2  and  3  will  be  briefly  reviewed. 
As  indicated  in  Figure  2,  a  cracked  body  subjected  to  arbitrary  elastic 
loading  (member  A)  can  be  resolved  into  two  components  (B  and  C) .  Member 
B  consists  of  the  uncracked  geometry  loaded  as  in  case  A.  A  stress  free 
line  is  introduced  along  the  desired  crack  plane  by  canceling  the  resulting 
hoop  stress  through  application  of  -p(x/R) .  Here  x  is  the  distance  measured 
from  the  edge  of  the  hole  (its  radius  is  R)  and  p(x/R)  is  the  hoop  stress 
distribution  along  the  crack  plane  caused  by  loading  the  unflawed  member. 

Since  the  crack  plane  is  now  stress  free,  a  crack  can  be  introduced.  Member 
C  consists  of  the  cracked  body  loaded  along  the  crack  faces  with  the  pressure 
+p(x/R).  By  superposition,  the  stress  intensity  factor  for  member  A  is  the 
linear  sum  of  the  stress  intensity  factors  for  members  B  and  C. 


ka  -  kb  +  Kc 


(1) 


Now,  since  the  crack  plane  in  member  B  is  stress  free.  Kg  =  0,  and  KA  “  Kc* 
Thus,  the  stress  intensity  factor  for  the  original  problem  is  identical  to 
that  for  a  crack  loaded  with  the  unflawed  hoop  stress  distribution. 

If  the  crack  face  pressure  is  defined  in  terms  of  a  polynomial  expansion, 
it  is  possible  to  obtain  a  "general"  solution  in  terms  of  the  polynomial 
coefficients.  This  latter  point  is  demonstrated  by  the  superposition  described 
in  Figure  3,  where 

p(x/R)  -  AQ  +  A1(x/R)  +  A2(x/R)2  +  ...  (2) 

Here  the  crack  face  pressure  is  expressed  in  terms  of  the  dimensionless 
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distance  x/R  from  the  edge  of  the  hole,  and  the  (n  *  0,  1,  2,  ...) 
are  the  polynomial  coefficients  determined  for  the  unflawed  hoop  stress 
distribution  (the  have  units  of  stress) .  Note  that  the  stress  intensity 
factors  for  the  cracks  shown  in  Figure  3  are  again  linearly  superimposed. 

Kj  -  KQ  +  Kx  +  K2  +  . . .  =  2Kn  (3) 

Here  is  the  mode-one  stress  intensity  factor  for  the  complete  crack  face 
loading  (equivalent  to  the  original  stress  intensity  factor  for  member  A 
in  Figure  2),  and  the  are  the  stress  intensity  factors  due  to  the  individual 
pressure  terms  An(x/R)n  in  the  expansion  for  p(x/R) . 

Note  that  if  the  are  obtained  for  arbitrary  polynomial  coefficients 
A^,  and  are  expressed  in  the  dimensionless  form  given  by  Eq.  4,  the  results 
are  readily  applied  to  other  plate  loadings. 

K  (4) 

M  =  — 
n  A  \ffiF 
n 


Here  the  magnification  factor  M  is  defined  for  each  individual  term  in 

n  - 

the  crack  face  pressure  expansion  and  c  is  the  crack  length  measured  in 
the  radial  direction.  Now,  given  the  A  an<^  crack  length  for  the  new  problem 


of  interest,  the  desired  stress  intensity  factor  K^.  is  found  from  the  general 


results  (the  M  )  by  performing  the  summation  of  Eq.  3. 

n  K 

Kt  -  £K  =Vttc"  >  - F= —  A 

I  n  <£_ _  A  Vue  n 

n 


(5) 


=  Vttc  A  *  K_ 
n  n  I 

Thus,  the  magnification  factors  obtained  for  the  crack  face  pressure 
loading  represent  general  solutions  which  may  be  used  to  compute 
stress  intensity  factors  for  any  loading  once  the  unflawed  hoop  stress  is 
found  and  represented  by  Eq.  2.  Note  that  although  the  hoop  stress  must 
still  be  determined  for  the  unf la wed  fastener  hole  problem  of  interest, 
the  fact  that  the  cracked  hole  does  not  have  to  be  directly  analyzed  greatly 
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simplifies  computational  effort.  It  is  anticipated  that  the  unflawed 
problem  will  be  analyzed  by  a  conventional  stress  analysis  method  (e.g. 
finite  elements). 

This  general  analysis  procedure  is  demonstrated  in  subsequent  sections 

by  several  example  calculations  with  through-the-thickness,  quarter- 

elliptical,  and  semielliptical  fastener  hole  cracks.  These  examples 

employ  results  given  here  for  the  dimensionless  crack  face  pressure 

stress  intensity  factor  M^.  Prior  to  discussing  those  results,  however, 

the  analysis  procedures  used  to  obtain  the  crack  face  pressure  solutions 

(the  K  and  subsequently  M  )  are  described, 
n  n 
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SECTION  II 


CRACK  FACE  PRESSURE  ANALYSIS  METHODS 


This  section  describes  the  analysis  procedures  used  to  obtain  the 
general  crack  face  pressure  solutions  for  the  two-dimensional  through- 
the-thickness  and  the  three-dimensional  corner-  and  surface-cracked  hole 
configurations.  A  weight  function  approach  was  used  for  the  through-the- 
thickness  cracks,  while  the  finite  element  -  alternating  method  was  em¬ 
ployed  for  the  quarter-  and  semielliptical  cracks. 

THROUGH-THE-THICKNESS  CRACKS 

The  through-cracked  hole  problem  is  shown  in  Figure  4.  An  earlier 
solution  (5)  developed  for  the  arbitrary  crack  face  pressure  loading  was 
used  to  compute  for  the  individual  pressure  term  A^Cx/R)11.  As  described 
in  Reference  5,  weight  function  methods  allowed  computation  of  the  crack 
face  pressure  solution  from  stress  intensity  factor  (6)  and  crack  opening 
displacement  results  for  remote  tt  sile  loading. 

The  crack  face  pressure  solution  obtained  by  the  weight  function  method 


is  given  by 


/  p(x/R) 


dx 


(6) 


Here  K*  is  the  stress  intensity  factor  for  remote  tensile  loading  applied 
perpendicular  to  the  crack  plane  (6),  c  is  the  crack  length,  R  is  the 
hole  radius,  T)  is  the  crack  opening  displacement  profile  corresponding 
to  K*,  p(x/R)  is  the  crack  face  pressure  as  before,  and  H  is  a  material 
constant  defined  in  terms  of  the  elastic  modulus  and  Poisson's  ratio. 


A  procedure  for  computing  the  derivative  of  the  crack  opening  profile 
with  respect  to  crack  length  (^J)  is  described  in  Reference  5.  A  com¬ 
puter  program  was  employed  to  numerically  integrate  Eq.  6  for  the  crack 
face  pressure  loading  p(x/R)  defined  by  Eq.  2.  Calculations  reported 
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in  References  4  and  5  indicate  that  the  stress  intensity  factor  solutions 

obtained  by  applying  Eq.  6  with  the  linear  superposition  procedure  outlined 

earlier  in  conjunction  with  Figure  2  give  excellent  stress  intensity  factor 

results  when  the  unflawed  hoop  stress  distribution  is  used  to  define  the 

crack  face  pressure  solution.  Thus,  based  on  this  successful  experience 

with  the  arbitrary  crack  face  pressure  solution,  Eq.  6  was  used  to  compute 

K  (and  the  dimensionless  magnification  factor  M  )  for  the  crack  face 
n  n 

loading  A  (x/R)n.  These  results  for  M  are  summarized  in  a  later  section 
n  n 

of  this  report. 

CORNER  AND  SURFACE  CRACKS 

The  quarter-elliptical  and  semielliptical  surface  cracks  located  at 
the  fastener  holes  are  shown  in  Figures  5  and  6.  Here  T  is  the  plate  thick¬ 
ness,  R  is  the  hole  radius,  a  is  the  semiaxis  of  the  elliptical  crack  bound¬ 
ary  measured  along  the  bore  of  the  hole,  c  is  the  semiaxis  length  in  the 
radial  direction,  and  $  is  the  "parametric"  angle  which  defines  points  along 
the  border  of  a  circle  circumscribed  around  the  part  elliptical  crack.  A 
projection  from  the  circle  defines  points  along  the  crack  border.  The  x  and 
z  coordinates  (see  Figures  5  and  6)  of  points  along  the  crack  perimeter  are 
given  by  x  =  c  cos  $  and  z  -  a  sin  $ . 

The  finite  element  -  alternating  method  was  used  to  obtain  the  stress 
intensity  factor  variation  around  the  perimeter  of  the  corner-  and  surface- 
cracked  holes.  The  finite  element  -  alternating  method  was  developed  by 
F.  W.  Smith  and  his  associates  and  has  been  successfully  applied  to  several 
surface  crack  problems  (7-10).  Specific  use  was  made  here  of  computer  codes 
developed  by  Smith  and  Kullgren  for  the  analyses  described  in  Reference  10. 

The  finite  element  -  alternating  method  Involves  iterative  superposition 
of  a  three-dimensional  finite  element  solution  for  an  uncracked  body  subjected 
to  a  prescribed  surface  loading  (11)  and  a  solution  for  a  flat  elliptical 
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crack  in  an  infinite  body  loaded  with  a  prescribed  nonuniform  surface 
pressure  (12) .  Iteration  between  these  two  solutions  approximates  the 
surface  crack  boundary  conditions  and  provides  mode-one  stress  intensity 
factors  and  crack  opening  displacements  for  the  three-dimensional  flaw 
geometry.  Additional  details  of  the  finite  element  -  alternating  method 
are  presented  in  References  7-10. 

Comparisons  of  stress  intensity  factors  obtained  by  the  finite  element  - 
alternating  procedure  with  other  experimental  and  numerical  solutions  for 
cracked  fastener  holes  indicate  that  the  analysis  method  provides  excellent 
stress  intensity  factor  results  for  this  class  of  problem.  Kullgren  and  Smith 
(8-10),  for  example,  have  applied  the  finite  element  -  alternating  method 
to  a  series  of  surface- cracked  holes  in  large  plates.  Results  are  presented 
for  three  crack  locations  near  fastener  holes  in  plates  with  two  hole  diameter 
to  plate  thickness  ratios.  Plate  loading  conditions  included  remote  uni¬ 
axial  tension  and  100  percent  transfer  of  load  to  a  fastener  filled  hole. 

Those  results  have  been  verified  in  several  ways.  Three-dimensional  calcula¬ 
tions  (8-10)  by  the  finite  element  -  alternating  method  for  through-cracked 
holes  loaded  in  remote  tension  agree  within  3.2  percent  for  single  crack 
results  reported  by  Tweed  and  Rooke  (13)  and  within  5  percent  of  Bowie's 
double  flaw  solution  (6) .  Corner  crack  results  agree  to  within  about  10 
percent  of  finite  element  data  reported  by  Raju  and  Newman  (14).  Other 
favorable  comparisons  of  the  finite  element  -  alternating  method  are  reported 
in  Reference  10  for  experimental  data  obtained  for  static  fracture  (15)  and 
fatigue  crack  growth  tests  (16)  with  corner- cracked  holes.  A  more  recent 
study  (17)  involving  simulated  aircraft  spectrum  loading  with  2124-T851  alu¬ 
minum  specimens  indicated  that  stress  intensity  factors  obtained  by  the  current 
analysis  accurately  predict  fatigue  crack  growth  life  and  changes  in  crack 
shape. 
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In  view  of  this  prior  success  with  the  finite  element  -  alternating 
method,  this  analysis  technique  was  used  to  determine  the  crack  face  pressure 
stress  intensity  factor  solutions  for  the  corner-  and  surface-cracked  holes. 
Those  results,  along  with  the  through-cracked  hole  data  obtained  by  the 
weight  function  method  are  discussed  in  the  following  section. 


SECTION  III 


STRESS  INTENSITY  FACTOR  RESULTS 


The  objective  of  this  section  is  to  present  the  stress  intensity  factor 
solutions  obtained  for  the  flawed  fastener  holes.  Application  of  these 
"general"  crack  face  pressure  results  to  several  example  problems  are 
discussed  in  the  following  section. 

THROUGH-CRACKED  HOLES 

As  described  earlier,  a  weight  function  method  was  used  to  compute  stress 

intensity  factors  for  through-cracked  holes  loaded  with  a  crack  face  pressure 

A^  (x/R)n.  Seven  individual  pressure  terms  were  considered  as  n  varied  in 

the  range  0<^n<_6.  Although  both  single  and  double  symmetric  cracks  were 

examined,  the  dimensionless  crack  lengths  were  restricted  to  the  range 

0.1<c/R<2.5.  The  stress  intensity  factor  results  are  summarized  in  the 

dimensionless  form  M  =  K  /(A  Vttc)  in  Tables  1  and  2.  These  data,  along 

n  n  n 

with  similar  results  for  through-cracked  rings,  were  reported  earlier  in 
Reference  18. 

CORNER -CRACKED  HOLES 

The  quarter-elliptical  corner  crack  geometry  is  shown  in  Figure  5. 

The  cracked  hole  problems  considered  here  are  summarized  in  Table  3.  In 
all  cases,  the  plate  hole  diameter  to  thickness  ratio  2R/T  was  fixed  at 
1.0,  Poisson's  ratio  was  kept  at  a  value  of  0.25,  and  the  "large"  circular 
plate  had  an  outer  radius  “  12  R. 

Although  solutions  for  a  large  matrix  of  crack  shapes  would  be  desirable, 
the  alternating  method  requires  fairly  large  computer  costs.  Thus,  care 
was  taken  to  select  the  problems  considered.  The  US AF  damage  tolerance 
specifications  described  in  Reference  1  require  analysis  of  quarter  circular 
corner  cracks  (a/c  =  1.0).  As  discussed  in  Reference  10,  however,  the 
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computer  codes  employed  by  the  present  authors  have  a  numerical  instability 
at  a/c  =  1.0.  Although  this  instability  is  associated  with  the  computer 
programming,  and  has  no  physical  significance,  it  prevented  direct  analysis 
of  quarter  circular  corner  cracks.  Thus,  aspect  ratios  a/c  *  0.9  and  1.1 
were  selected  to  closely  bracket  the  crack  shapes  specified  in  the  USAF 
requirements.  Emphasis  was  placed  on  results  for  a/c  »  1.1  rather  than 
0.9  because  experience  indicates  that  naturally  occurring  fatigue  cracks 
often  have  a/c  values  between  1.0  and  1.5  (16). 

The  loading  conditions  considered  here  included  the  first  four  individual 
terms  in  the  polynomial  expansion  for  the  crack  face  pressure.  Thus,  n 
covered  the  range  0«ns3  for  the  pressure  distribution  given  by  A^Cx/R)11. 
Computations  were  also  performed  for  a  remote  tensile  stress  a  applied 
perpendicular  to  the  crack  plane  (see  Figure  7) . 

The  corner-cracked  hole  stress  intensity  factor  results  are  summarized 
in  Tables  4  and  5.  Here  the  dimensionless  magnification  factor  is  presented 
as  a  function  of  crack  shape  (a/c),  size  (a/t) ,  position  along  the  crack 
front  <P  ,  and  applied  load.  These  results  were  presented  earlier  in  Reference 
19. 


The  magnification  factor  in  Table  4  is  defined  by 


M 


A  \/jfa 


(7) 


n 

As  described  earlier,  is  the  stress  intensity  factor  caused  by  An(x/R)n 
in  the  crack  face  pressure  expansion  defined  by  Eq.  2.  Note  that  since 
the  polynomial  coefficients  A^  have  units  of  stress  (see  Eq.  2),  the  magni¬ 
fication  factor  M  is  dimensionless.  The  data  in  Table  4  are  for  values 
n 

of  n  ■»  0,  1,  2,  and  3. 

The  results  in  Table  5  are  for  a  remote  stress  <r  applied  perpendicular 
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where  is  the  stress  intensity  factor  caused  by  the  remote  stress  <r. 

All  of  the  results  presented  in  Tables  4  and  5  are  for  a  material  having 
a  Poisson's  ratio  of  0.25,  a  hole  diameter  to  plate  thickness  ratio  2R/T  » 
1.0  and  are  reported  for  30  degree  increments  of  the  "parametric"  angle  <p 
defined  in  Figure  ■*. 

In  contrast  to  the  crack  face  pressure  results  (Table  4),  the  remote 
stress  (Table  5)  data  include  through-the-thickness  effects  in  the  stress 


intensity  factor  calculations.  Reference  20  discusses  a  three-dimensional 


stress  analysis  for  remotely  loaded  holes  in  plates  of  arbitrary  thickness, 
and  indicates  that  the  unflawed  hoop  stress  for  the  geometry  considered 
here  (hole  diameter  =  plate  thickness)  can  vary  by  as  much  as  10  percent 
in  the  thickness  direction  (z  direction  in  Figures  5  and  6)  at  the  bore  of 
the  hole.  As  described  in  Reference  10,  direct  application  of  the  finite 
element  -  alternating  method  to  the  remote  load  problem  (the  Table  5  results) 
takes  this  three-dimensional  effect  into  account. 


SURFACE -CRACKED  HOLES 

The  semielliptical  surface~cracked  hole  geometry  is  shown  in  Figure  6. 
The  flaw  shapes  and  loading  conditions  considered  here  are  summarized  in 
Table  6.  The  loading  conditions  again  consisted  of  the  individual  crack 
face  pressure  terms  A^(x/R)n  where  0<ns3,  and  a  remote  tensile  stress 
a  applied  perpendicular  to  the  crack  plane  as  in  Figure  7.  Crack  shapes 
include  aspect  ratios  a/c  =  1.11,  1.5,  and  2.0  and  crack  lengths  2a/T  * 
0.2,  0.4,  0.6,  0.8,  and  1.0. 

The  crack  shapes  studied  were  selected  to  agree  with  those  found  to 
occur  "naturally"  for  this  flaw  configuration.  Experiments  with  surface- 


cracked  hole  specimens  loaded  in  remote  tension  reported  by  Kullgren  and 
Smith  (15),  for  example,  gave  crack  shapes  in  the  range  of  2.17fa/c<1.52 
with  a/c  =  1.85  as  the  average  value.  Fractographic  examination  of  cracked 
fastener  holes  in  retired  jet  engine  turbine  disks  reported  in  Reference  21 
give  typical  aspect  ratios  on  the  order  of  a/c  =  1.4.  Specimens  prepared  for 
a  nondestructive  inspection  study  (22)  associated  with  a  recent  damage  toler¬ 
ance  analysis  conducted  on  another  USAF  turbine  engine  gave  small  semi¬ 
elliptical  cracks  with  aspect  ratios  in  the  range  of  1.2fa/c<1.0.  Hsu  et.  al. 
(23)  also  report  fatigue  'at a  for  fastener  holes  with  two  symmetric  surface 
cracks  (one  on  each  .vide  of  the  hole)  with  aspect  ratios  1.50<a/c<1.25. 

Dimensionless  stress  intensity  factors  for  the  two  loading  configurations 
(crack  face  pressure  and  remote  stress)  are  presented  in  Tables  7  and  8  as  a 
function  of  position  along  the  crack  perimeter  (<j>) ,  crack  size  (2a/T) ,  and 
flaw  shape  (a/c).  Again  the  plate  hole  to  thickness  ratio  2R/T  =  1.0, 

Poisson's  ratio  has  a  value  of  0.25,  and  the  remote  stress  loading  results 
incorporate  the  same  through-the-thickness  variation  in  stress  described 
earlier  for  the  corner-cracked  hole  configuration.  The  magnification  factors 
M  and  M  are  defined  by 


M  = 


n 


K 

_ n _ 

An/iTc 


^9) 


and 


M  =  ~~  (10) 

a/irc 

These  results  are  also  reported  in  Reference  24.  Note  that  the  magnification 
factor  for  the  semielliptical  surface  cracks  (Eqs.  9  and  10)  are  defined  with 
the  dimension  "c"  in  the  denominator  rather  than  "a"  as  used  in  Eqs.  7  and  8 
for  corner  cracks. 
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STRESS  FUNCTION  COEFFICIENTS 

Tables  4,  5,  7  and  8  present  tabulated  stress  intensity  factors  for  the 
corner  and  embedded  (surface)  cracks  at  30°  increments  of  the  parametric  angle  <p. 
The  original  computer  output  gave  results  for  other  angles  and  also  provided 
crack  opening  displacement  information.  Although  a  complete  tabulation  of  all 
the  available  data  is  too  voluminous  for  this  report,  it  is  possible  to  provide 
information  which  can  be  used  to  reconstruct  both  stress  intensity  factors  and 
crack  opening  displacements  for  additional  values  of  the  parametric  angle  <f>. 
Appendix  A,  for  example,  summarizes  a  set  of  10  stress  function  coefficients 
which  were  computed  (see  Refs.,  10  and  12  for  details)  for  each  crack  configuration 
considered  here.  Computer  programs  described  in  Appendices  B  and  C,  use 
these  stress  function  coefficients  to  recalculate  stress  intensity  factors 
and  crack  opening  disp lacemen ts  as  a  function  of  the  parametric  angle  <f>.  Thus, 
as  demonstrated  in  the  appendices,  the  reader  can  readily  obtain  complete 
descriptions  of  stress  intensity  factors  and  crack  opening  displacements  for 
all  of  the  problems  summarized  in  Tables  3  and  6.  (These  problems  are  relisted 
and  numbered  in  Tables  A-l  and  A-2  of  Appendix  A.) 
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SECTION  IV 


EXAMPLE  PROBLEMS 


The  purpose  of  this  section  is  to  demonstrate  application  of  the  "general" 
stress  intensity  factor  solutions  reported  in  Tables  1,  2,  4,  and  7.  These 
crack  face  pressure  results  are  applied  to  several  example  problems  in  order 
to  further  explain  the  calculations  of  Eq.  5,  and  to  indicate  the  general 
utility  of  the  tabulated  data.  Where  possible,  these  calculations  are  com¬ 
pared  with  independent  results  in  order  to  assess  the  accuracy  of  the  super¬ 
position  procedure. 


EXAMPLE  1;  THROUGH-CRACKED  HOLE  LOADED  IN  REMOTE  TENSION 

This  example  demonstrates  the  use  of  crack  face  pressure  results  in 
Table  1  to  compute  the  stress  intensity  factor  for  a  through-cracked  hole 
loaded  in  remote  tension  perpendicular  to  the  crack  plane.  This  flaw  geometry 
was  analyzed  by  Bowie  (6)  and  stress  intensity  factor  results  are  reported 
in  Reference  25.  For  demonstration  purposes,  assume  crack  length  c  =  0.5R 
and  the  applied  stress  is  <r. 

The  first  step  in  the  computational  procedure  is  to  determine  the  unflawed 

hoop  stress  v  caused  by  the  remote  load.  A  two-dimensional  stress  analysis 
rl 

for  <r„  is  given  in  Reference  26.  Least  squares  procedures  were  employed 
n 

here  to  express  that  result  in  the  polynomial  expansion  given  by  Eq.  2. 

The  following  expression  agrees  within  three  percent  of  the  exact  analysis 
for  0<x/Rsl.0,  where  <r  is  the  remote  stress. 

crH  =<jj^2.9190  -  5.2853(x/R)  +  6.2574(x/R)2  -  2 . 6910(x/R) 3 j  (11) 

Now,  letting  the  crack  face  pressure  p(x/R)  «  cr  ,  the  pressure  polynomial 
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coefficients  (A  )  and  magnification  factors  (M  )  from  Table  1  (for  c/R  =  0.5) 


Ao 

=  2.9190 

Mo  = 

.9119 

A1 

=  -5.2853 

a* 

Mi 

.3020 

A2 

=  6.2574 

M2  = 

.1207 

A3 

=  -2.6910 

<r 

N 

.05167 

The  desired  stress  intensity  factor  is  now  given  by  Eq.  5  as 
K  =  t  K  =  VttF)  -  -n-  n-  =  Va  M 

1  n  L  n  n 

n  =  0  n 

=  «r  ^/5rc  [(2.9190)  (.9119)  -  (5.2853)  (.3020)  +  (6.2574)  ( .1207)  - 

(2 . 6910)  ( .  05167)]  =  1.6819  <rv^c  =  ^ . 


The  corresponding  result  given  in  Reference  25  (for  c/R  ■  0.5)  is 
K  -  1.73  Thus,  the  solution  obtained  by  the  superposition  procedure 

agrees  within  three  percent  of  the  "exact"  result  found  by  direct  analysis 
(6)  of  the  cracked  member. 


EXAMPLE  2:  THROUGH-CRACKED  COLDWORKED  HOLE 

Coldworking  fastener  holes  is  an  effective  method  for  extending  the 
fatigue  lives  of  mechanically  fastened  joints.  The  process  involves  intro¬ 
ducing  residual  compressive  hoop  stresses  at  the  hole  edge  by  pulling  an 
oversize  mandrel  through  the  hole.  These  residual  stresses  reduce  the 
effective  stress  concentration  when  remote  load  is  applied,  and  can  greatly 
extend  the  fatigue  life  of  fastener  holes  (2-4). 

The  purpose  of  this  example  is  to  demonstrate  application  of  the  tabu¬ 
lated  crack  face  pressure  solutions  to  the  coldworking  process.  As  discussed 
in  Reference  4>  consider  a  0.261  inch  (6.63  mm)  diameter  hole  in  a  0.261  inch 
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(6.63  mm)  thick  plate  of  7075-T651  aluminum.  A  0.274  inch  (6.96  mm)  diameter 
mandrel  is  pulled  through  the  fastener  hole  and  causes  a  residual  stress  field 
around  the  hole.  The  plate  is  then  loaded  with  a  remote  uniaxial  stress  of 
40  ksi  (276  MPa).  The  objective  here  is  to  compute  for  a  through-the- 

thickness  crack  length  c/R  =  0.8  (c  =  .104  inch  =  2.65  mm). 

In  order  to  apply  the  superposition  calculation  to  the  coldworked  hole 
configuration,  it  is  again  necessary  to  determine  the  hoop  stress  distribu¬ 
tion  for  the  uncracked  case.  One  may,  for  example,  establish  the  coldworked 
hole  residual  stress  field  by  the  "exact"  elastic-plastic  solution  given  by 
Hsu  and  Forman  (27)  for  an  infinite  sheet  containing  a  pressurized  circular 
hole.  Their  plane  stress  result  is  based  on  the  deformation  theory  and 
modified  Ramberg-Osgood  law.  A  small  computer  program  was  written  here  to 
assist  in  solving  the  fairly  involved  equations  which  represent  the  Hsu-Forman 
solution  and  used  to  obtain  the  residual  stress  distributions  shown  in  Figure 
8-  (Note  that  there  are  several  typographical  errors  in  the  equations  given 
in  Reference  27). 

Figure  8  indicates  that  large  compressive  hoop  stresses  are  left  next  to 
the  hole  edge  by  the  pre-expansion  and  are  balanced  further  away  from  the  hole 
by  a  corresponding  tensile  field.  Since  large  tensile  stresses  are  concen¬ 
trated  at  the  edge  of  an  unworked  hole  when  the  plate  is  loaded  in  uniaxial 
tension,  the  residual  stress  distribution  reduces  the  peak  tension  next  to 
the  hole.  An  opposite  effect,  which  could  lead  to  more  plastic  flow,  and  to 
a  change  in  the  residual  stress  pattern,  would  occur  if  the  plate  were  loaded 
in  compression  following  coldworking. 

For  example,  the  residual  hoop  stress  in  Figure  8  was  fit  to  within  three 
percent  by  a  polynomial  as  before,  and  gave  the  following  result. 
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n  =  0 

=-130.53  +  296.25(x/R)  -  226.67(x/R)2  +  72.936(x/R)3 
Here  the  polynomial  coefficients  have  units  of  ksi  and  x/R  is  restricted 
to  the  range  0  <x/Ri  1.0.  Now,  the  unflawed  hoop  stress  for  the  remotely 
loaded  coldworked  hole  is  found  by  adding  Eqs.  Hand  12.  Letting  o  =  40  ksi 
in  Eq.  11  (the  remote  tensile  stress  of  interest)  gives 
o’  =  o’  -f*  cr 

total  remote  coldwork  (13) 


=  S  (0-A  +B)(x/R)n=  V  C  (x/R) n 
Z-r  n  n  /.  n 


n  =  0 


n  =  0 


=-13.770  +  84.838(x/R)  +  23.626(x/R)2  -  34.704(x/R)3 
Note  that  the  superposition  of  Eq.  13  implies  that  the  remote  stress  field 
is  applied  elastically  and  does  not  alter  the  residual  stresses  introduced 
by  coldworking.  Since  the  sum  of  the  remote  and  residual  stresses  never 
exceeds  yield,  this  superposition  is  valid  for  the  present  example. 

Now,  the  desired  stress  intensity  factor  for  the  0.104  inch  j.ong  through- 
crack  (c/R  =  0.8)  can  be  computed  by  superposition  as  before.  Table  1  and 
Eq.  13  provide  the  following  information. 


MQ  =  .8610 
Kj,  =  .4693 
M2  =  .3030 
M3  =  .2089 

The  summation  of  Eq.  5  gives 


CQ  =  -13.770  ksi 
C1  =  84.838  ksi 
C2  =  23.626  ksi 
C3  =  -34.704  ksi 


Kt  K  =  \Zic  Vh  C 
I  L  n  /.  n  n 

*  v/Or)(.1044)  {(.8610) (-13.770)  +  ( .4693) (84 .838)  +  ( .3030) (23 .626)  + 
(.2089) (-34.704)}=  15.96  ksi  -  in  1/2  =  K^. 
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This  result  is  plotted  as  the  large  "star"  on  Figure  9,  and  is  compared 
with  experimental  measurements  on  coldvorked  hole  specimens  reproduced  from 
Reference  4.  The  line  labeled  "Coldwork"  was  obtained  by  superposition 
calculations  similar  to  that  employed  here  except  that  a  higher  order 
polynomial  fit  was  used  to  extend  the  x/R  range  in  Eqs.  11-13.  The  "zero  R" 
values  indicated  for  small  crack  lengths  actually  resulted  from  calculations 
which  gave  negative  stress  intensity  factors  (note  that  Eq.  13  gives  nega¬ 
tive  hoop  stresses  for  small  x/R).  Since  negative  K's  have  no  physical 
meaning,  all  negative  results  were  set  equal  to  zero  in  Fig.  9.  The  line 
labeled  "Bowie"  is  for  the  unworked  case  (remote  load  only) .  The  experi¬ 
mental  data  were  obtained  from  fatigue  crack  growth  rate  measurements  in 
precracked  coldworked  hole  specimens  (see  Reference  4  for  additional  details). 
Note  that  the  predicted  stress  intensity  factors  agree  well  with  the  experi¬ 
mental  data,  indicating  the  superposition  procedure  can  give  good  predic¬ 
tions  for  the  coldworking  process.  Similar  results  are  given  in  Reference 
4  for  other  remote  load  applications,  and  indicate  that  the  stress  intensity 
factor  calculations  can  provide  good  estimates  of  the  fatigue  growth 

lives  of  coldworked  holes. 

EXAMPLE  3;  CORNER -CRACKED  HOLE,  REMOTE  TENSION 

Consider  a  corner  crack  emanating  from  an  open  hole  as  in  Example  1. 
Stress  intensity  factors  for  this  flaw  configuration  were  obtained  directly 
by  the  finite  element  -  alternating  method  and  are  reported  in  Table  5. 

The  goal  of  this  example  is  to  demonstrate  application  of  the  crack  face 
pressure  results  in  Table  4  by  employing  them  to  Independently  calculate 
the  remote  load  results  (Table  5)  through  the  superposition  procedure. 

The  specific  objective  is  to  compute  Kj/(o/rta)  at  4>  =*  90°  for  a  crack  which 
has  an  aspect  ratio  a/c  -1.5  and  length  a/t  *  0.8. 

The  unflawed  hoop  stress  distribution  is  again  given  by  Eq.  11,  and 
Table  4  provides  the  appropriate  crack  face  stress  intensity  coefficients  Mq. 
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Mq  =  .4994 

Ao 

=  2.9190  <t 

*  .0939 

A1 

=  -5.2853  tr 

>i2  =  .0280 

A2 

=  6.2574a- 

M3  =  .0027 

A3 

=  -2.6910  a 

Equation  5  now  gives 

Kr  =  <rV£a  {(.4994)  (2.9190)  +  (.  0939)  (-5. 2853)  +  ( .0280)  (6.2574)  + 

(.0027)  (-2. 6910)}  =  1.1294  oVttT 

The  corresponding  result  obtained  by  an  independent  application  of  the 
finite  element  -  alternating  method  as  described  in  Reference  10  and  reported 
in  Table  5,  is  K^(<r>/?a)  =  1.1660  (a  3.1  percent  difference).  This  difference 
could  be  due  to  the  error  in  fitting  the  exact  hoop  stress  solution  by  Eq.  11 
or,  as  discussed  earlier,  to  through-the-thickness  effects,  or  to  approxi¬ 
mations  employed  by  polynomial  fitting  routines  in  the  finite  element  - 
alternating  method. 


EXAMPLE  4:  SURFACE -CRACKED  HOLE,  REMOTE  TENSION 

Consider  another  plate  with  an  open  hole  loaded  in  remote  tension.  For 
this  example,  assume  a  surface  crack  occurs  along  the  bore  of  the  hole  (see 
Figure  6),  and  compute  the  variation  in  stress  intensity  factor  around  the 
perimeter  of  a  crack  which  has  an  aspect  ratio  a/c  =  1.11  and  length  2a/T  « 
0.8.  The  unflawed  hoop  stress  is  again  given  by  Eq.  11,  and  the  stress 
intensity  factors  are  now  found  in  Table  7.  At  the  crack  perimeter  location 
defined  by  =  90°  (where  the  crack  intersects  the  bore  of  the  hole),  the 


M  are 
n 


MQ  -  0.6836 
Mx  -  0.1030 
M2  -  0.0270 
M3  -  0.0061 
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Computing  the  stress  intensity  factor  by  Eq.  5  gives 

Kt  =  JOc  V  AM 

I  n  /  n  n 

n  =  0 

=  ar\firc  {(2.9190)(0.6836)  -  (5.2853)  (0. 1030)  +  (6.2574)  (0.0270)  - 
(2.6910)(0.0061)}  =  1 . 6036  <r  S^Tc 

Note  that  the  result  obtained  directly  by  the  finite  element  -  alternating 
method  reported  in  Table  8  for  this  case  (a/c  =  1.11,  2a/T  =  0.8,  <p  =  90°) 
is  Kj.  =  1.625  ary/ffc.  Similar  calculations  for  <P  =  0°  ,  30°,  and  60°  yield 
respective  values  of  0.995,  1.113,  and  1.401  for  K^ir\/Sc.  These  results 
also  compare  quite  favorably  with  the  independent  calculations  reported 
in  Table  8. 

For  comparison  purposes,  the  superposition  procedure  described  in 
Examples  3  and  4  was  used  to  recompute  all  corner- cracked  and  surfaced- 
cracked  hole  results  reported  in  Tables  5  and  8.  In  this  manner,  it  was 
possible  to  determine  the  accuracy  achieved  by  the  superposition  method 
for  surface  and  corner  cracks.  For  the  surface  cracked  holes,  the  Table 
8  results  (remote  loading  applied  directly)  agreed  within  5.3  percent  of 
all  stress  intensity  factors  found  from  the  c rack  face  pressure  solutions. 
The  largest  discrepancy  occured  at  4>  =  90°  for  the  smallest  cracks 
(2a/T  =  0.2)  for  all  three  flaw  shapes  (a/c  =  1.11,  1.5,  and  2.0).  For 
the  corner- cracked  holes,  the  Table  5  results  agree  within  9  percent  of 
the  superposition  calculations  (see  Example  3)  for  all  cases.  The  largest 
difference  again  occurred  for  the  smallest  crack  lengths  (2a/T  *  0.2)  at  the 
point  where  the  crack  intersects  the  bore  of  the  hole  (<}>  =  0  in  Figure  5). 

In  general,  most  calculations  were  well  within  5  percent  agreement. 

Thus,  these  "calibration"  computations  indicate  that  the  superposition 
procedure  provides  an  effective  means  for  estimating  the  stress  intensity 
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factor  distribution  around  the  perimeter  of  corner-  and  surface-cracked 
holes.  The  fact  that  the  greatest  "error"  occurs  for  the  smallest  crack, 
lengths  along  the  bore  of  the  hole  is  consistent  with  the  thickness  varia¬ 
tion  in  hoop  stress  (20)  described  earlier. 


EXAMPLE  5:  CORNER -CRACKED  COLDWORKED  HOLES 


For  the  final  example,  consider  a  corner  crack  which  occurs  at  the 
coldworked  hole  described  in  Example  2.  Assume  the  crack  shape  a/c  =  1.1, 
the  flaw  length  a/T  =  0.4  (a  =  0.104  in  =  2.65  mm),  and  the  remote  applied 
stress  is  45  ksi  (310  MPa) .  The  uncracked  hoop  stress  can  again  be  found 
by  adding  the  remote  and  residual  coldworked  hole  stresses  as  in  Eq.  13 
(in  this  case,  <r  =  45  ksi).  For  this  example,  the  Eq.  13  calculation  and 


Table  4  gives  the  following  for  4>  =  0. 


Mq  =  0.6656 
Mx  =  0.3329 
M_  *  0.1943 
M3  =  0.1202 

The  superposition  calculation  now  gives 
3  3 


CQ  =  0.8250  ksi 
C1  =  58.41  ksi 
C2  =  54.91  ksi 
C3  =  -48.16  ksi 


V  K  =  y/na  S  M  C 

L  n  L  n  n 


=  14.22  Ksi-in 


n  =  0 


n  =  0 


Similar  calculations  at  other  points  around  the  crack  border  (<f>  =  30°,  60°, 

and  90°)  give  K^  values  of  12.69,  7.15,  and  3.16  Ksi-in^^ , 

Note  that  for  an  identical  crack  in  a  plate  which  is  not  coldworked, 

but  subjected  to  the  remote  tensile  stress  only,  the  results  of  Table  5 

(or  approximate  calculations  with  the  Table  4  crack  face  pressure  solutions 

1/2 

as  in  Example  5)  give  values  of  26.4,  26.0,  31.3  and  37.2  Ksi-in 

for  the  respective  parametric  angles  $  =  0°,  30°,  60°,  and  90°.  The  large 


reduction  in  K^  for  the  coldworked  hole  case  is  consistent  with  the  experimental 


\ 
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fact  that  coldworking  can  significantly  reduce  fatigue  crack  growth  rates 

at  fastener  holes.  Note  that  the  reduction  in  is  greatest  at  the  bore 

of  the  hole  < <j>  »  90°)  where  coldworking  lowers  to  3.16  Ksi-in^^  from 

1/2 

37.2  Ksi-in  ,  indicating  that  coldworking  is  most  effective  next  to  the 
hole  edge. 
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SECTION  V 


SUMMARY  AND  CONCLUSIONS 

The  objective  of  this  report  is  to  present  a  set  of  general  stress 
intensity  factor  solutions  which  may  be  used  to  analyze  practical  fastener 
hole  problems.  Results  are  given  for  radial  through- the- thickness,  quarter- 
elliptical  corner,  and  semielliptical  surface  (embedded)  cracks  which  are 
loaded  with  crack  face  pressure.  Through  linear  superposition,  these 
tabulated  results  are  readily  applied  to  other  types  of  loading.  This 
superposition  procedure  is  demonstrated  through  several  sets  of  example 
calculations. 

It  is  felt  that  the  dimensionless  stress  intensity  factors  reported 
in  Tables  1,  2,  4,  and  7  provide  an  effective  means  for  estimating  K^.  for 
practical  fastener  hole  crack  configurations.  In  applying  these  crack 
face  pressure  solutions  to  new  problems,  however,  the  designer  or  analyst 
should  consider  the  following  points. 

1.  The  crack  face  pressure  (unflawed  crack  plane  hoop  stresses)  should 
be  tension  over  the  entire  crack. 

2.  The  unflawed  hoop  stresses  must  be  adequately  described  by  the 
polynomial  expression  given  by  Eq.  2.  Through-the-thickness 
problems  can  employ  a  sixth  degree  fit,  while  the  corner  and 
surface  crack  geometries  are  limited  to  third  order  polynomials. 

Note,  however,  that  it  is  only  necessary  to  represent  the  hoop 
stress  over  the  range  0  X/R  <_  c. 
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3.  Through- the- thickness  variations  in  the  crack  face  pressure  are 
not  considered.  Thus,  caution  should  be  exercised  in  applying  the 
crack  face  pressure  solutions  to  situations  where  the  unflawed  hoop 
stress  exhibits  a  significant  through-the-thickness  variation.  The 
remote  stress  results  reported  in  Tables  4  and  8  do,  however,  contain 
through-the-thickness  effects. 

4.  The  corner  and  surface  crack  results  are  restricted  to  plates  which 
have  hole  diameters  equal  to  the  plate  thickness  (2R=T)  and  to 
materials  with  a  Poisson's  ratio  of  0.25. 

Questions  will  naturally  arise  regarding  extrapolation  of  the  present 

results  to  cases  of  differing  hole  diameters  to  Dlate  thickness  ratios, 
loading  conditions  and  crack  shapes  or  sizes.  Kullgren  and  Smith  (9.  10) 
list  four  effects  governing  the  magnitude  and  variation  of  mode-one 
stress  intensity  along  the  crack  border. 

1.  Hole  Stress  Concentrations.  Locations  on  the  crack  border  close 

to  the  hole  and  thus  in  a  region  of  elevated  stress  have  larger 
stress  intensity  factors  than  do  locations  more  removed. 

2.  Front  and  Back  Surfaces.  Crack  points  close  to  a  plate  front 

or  back  surface  or  near  a  crack  intersection  with  these  surfaces 
have  elevated  stress  intensity  factors. 

3.  Hole  Surface.  For  small  cracks  (relative  to  hole  diameter) 
the  hole  surface  has  an  effect  similar  to  crack  intersection 
with  the  plate  front  or  back  surfaces.  For  large  cracks, 
the  hole  restricts  crack  opening  resulting  in  lower  stress 

intensity  near  the  hole/crack  intersection. 

4.  Crack  Shape.  Stress  intensity  factors  are  highest  at  the  ends 

of  an  elliptical  crack  under  constant  load  (9). 
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In  general,  it  is  not  possible  from  the  present  results  to  separate 
or  determine  the  relative  importance  of  these  effects.  Any  change  to 
the  crack  or  plate  geometries  results  in  more  than  one  effect  operating. 
For  example,  earlier  work  (7,  9,  10)  indicates  a  reduction  of  2R/T 
from  1.0  to  0.5  while  all  other  parameters  remain  unchanged  resulted 
in  a  lowering  of  magnification  factors  at  all  locations  on  the  crack 
border  for  cases  of  remote  loading  (open  fastener  hole).  However, 
when  this  same  load  is  transferred  to  a  fastener  in  the  hole,  magnifi¬ 
cation  factors  increased  for  the  same  reduction  in  2R/T.  While  this 
surely  underscores  the  importance  of  the  Hole  Stress  Concentration 
Effect,  it  is  by  no  means  clear  that  for  other  plate/crack  geometries 
the  same  trend  will  be  seen.  In  the  absence  of  an  extension  of  the 
present  work,  the  authors  do  not  recommend  the  quantitative  extrapo¬ 
lation  of  the  results  presented  here. 
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Table  1.  Dimensionless  stress  intensity  factors  for  a  single-cracked  hole  loaded  with  a  crack 
face  pressure  p(x)  =  A  (x/R)n. 


Table  3  Summary  of  Corner- Cracked  Hole  Solutions 
(x  denotes  geometries  analyzed) 


Crack  Shape 

a/c 

Crack  Size 

a/T 

Loading 

Crack 

face 

pressure 

An(x/R)n 

Remote 

n=0 

n*=l 

n=2 

n=3 

0 

.5 

.2 

X 

X 

X 

X 

X 

.4 

X 

X 

X 

X 

X 

.5 

X 

X 

X 

X 

X 

.9 

.2 
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X 
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X 

X 

X 

1.1 

.2 

X 

X 

X 

X 

X 

.4 

X 

X 

X 

X 

X 

.6 

X 

X 

X 

X 

X 

.8 

X 

X 

X 

X 

X 

1.0 
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1.5 

.2 
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.4 
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X 

X 

X 

.6 
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X 

X 

X 

X 

.8 
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X 

X 

X 

X 

1.0 

X 

X 

X 

X 
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Table  4  (Cont'd)  Screes  Intensity  Magnification  Factors 
CM  ■  K  /A  /  sa)  for  Corner  bracked 
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Table  5  Stress  Intensity  Magnification  Factors 
for  Corner  Crack  Holes  Loaded  with  a 
Remote  Uniaxial  Stress  o 


a/c 

a/t 

M=K  /(a/lla) 

<J>=0 

o 

o 

CO 

60° 

O 

o 

.5 

.2 

.930 

1.064 

1.514 

1.890 

.4 

.954 

.931 

1.269 

1.762 

.5 

.982 

.929 

1.217 

1.696 

.9 

.2 

1.106 

1.220 

1.482 

1.640 

.4 

1.050 

.999 

1.235 

1.543 

.6 

1.060 

.944 

1.114 

1.421 

.8 

1.068 

.906 

1.041 

1.369 

1.1 

.2 

1.235 

1.312 

1.456 

1.539 

.4 

1.027 

1.009 

1.215 

1.447 

.6 

1.059 

.952 

1.090 

1.338 

.8 

1.038 

.896 

1.019 

1.292 

1.0 

1.121 

.899 

1.311 

1.947 

1.5 

.2 

1.263 

1.347 

1.383 

1.365 

.4 

1.070 

1.069 

1.190 

1.302 

.6 

1.050 

.969 

1.059 

1.210 

.8 

1.009 

.903 

.991 

1.166 

1.0 

1.465 

1.242 

1.750 

2.386 
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Table  7  Stress  Intensity  Magnification  Factors  for 
Semielliptical  Surface  Cracks  at  Holes  in  a 
Plate  Loaded  with  Crack  Face  Pressure  A  (x/R)n 
(2R  =  T  for  all  cases).  n 


a/c  2a/T 

n 

",  ■  V(\ 

/^) 

$  =  0° 

+  30° 

+  60° 

+  90° 

1.11  0.2 

0 

.6761 

.6800 

.6674 

.5973 

1 

.0882 

.0746 

.0448 

.0187 

2 

.0129 

.0098 

.0042 

.0011 

3 

.0020 

.0014 

.0004 

.0001 

0.4 

0 

.6336 

.6526 

.6697 

.6189 

1 

.1718 

.1467 

.0914 

.0427 

2 

.0507 

.0386 

.0168 

.0053 

3 

.0156 

.0111 

.0035 

.0005 

0.6 

0 

.6276 

.6465 

.6710 

.6361 

1 

.2561 

.2186 

.1375 

.0682 

2 

.1139 

.0868 

.0381 

.0128 

3 

.0528 

.0375 

.0118 

.0019 

0.8 

0 

.6453 

.6650 

.7003 

.6836 

1 

.3456 

.2957 

.1904 

.1030 

2 

.2040 

.1560 

.0706 

.0270 

3 

.1260 

.0898 

.0293 

.0061 

1.0 

0 

.6565 

.7250 

.8943 

1.0039 

1 

.4379 

.3863 

.2840 

.2022 

2 

.3222 

.2511 

.1278 

.0690 

3 

.2484 

.1791 

.0640 

.0215 
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Table  7  (Continued) 


Table  8  Stress  Intensity  Magnification  Factors  for 
Semi  elliptical  Surface  Cracks  at  Holes  in 
Large  Plates  Loaded  with  Uniaxial  Tensile 
Stress  a  (2R  =  T  in  all  cases). 


a/c  2a/T  M  =  Kj/faSZc) 


<j>  =  0° 

+30° 

o 

0 

o 

o 

l.ll 

0.2 

1.608 

1.687 

1.803 

1.733 

0.4 

1.244 

1.375 

1.619 

1.681 

0.6 

1.081 

1.204 

1.474 

1.620 

0.8 

1.021 

1.132 

1.413 

1.625 

1.0 

0.971 

1.166 

1.710 

2.206 

1.50 

0.2 

1.956 

1.970 

1.904 

1.702 

0.4 

1.630 

1.711 

1.786 

1.696 

0.6 

1.382 

1.491 

1.650 

1.654 

0.8 

1.275 

1.378 

1.564 

1.640 

1.0 

1.211 

1.403 

1.817 

2.095 

2.0 

0.2 

2.227 

2.192 

1.968 

1.620 

0.4 

1.994 

2.019 

1.913 

1.644 

0.6 

1.725 

1.801 

1.810 

1.632 

0.8 

1.554 

1.648 

1.723 

1.619 

1.0 

1.445 

1.629 

1.904 

1.944 
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Figure  1  Flawed  Fastener  Hole  Configurations 
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Figure  6:  Semielliptical  surface  (embedded)  crack 
located  along  bore  of  a  hole  in  a  large 
plate. 


<r 


cr 


re  7:  Cracked  hole  in  a  large  plate  loaded  with  remote  tensil 
stress  applied  perpendicular  to  the  crack  plane 


Residual  Stresses  (KSI) 


Distance  X  From  Edge  Of  Hole  (MM) 


X  ( Inches) 


Figure  8:  Residual  stress  distribution  due  to  coldworking  hole  as 
predicted  by  Hsu-Forman  solution 
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APPENDIX  A 


TABULAR  DATA 

1 .  INTRODUCTION 

Appendices  A  through  C  are  designed  to  provide  completeness  to  this 
report.  In  Appendix  A,  results  from  finite  element-alternating  method 
calculations  are  listed  for  all  150  cases  considered.  Appendix  B  de¬ 
scribes  two  short  computer  programs  which  can  be  used  to  regenerate  more 
comprehensive  descriptions  of  mode  one  stress  intensity  variations  along 
the  crack  border.  Appendix  C  contains  another  short  program  which  gives 
crack  opening  displacements.  The  computer  programs  of  Appendices  B  and 
C  use  the  tabular  results  of  Appendix  A  as  input . 

2.  DATA  DESCRIPTION 

Table  A-l  lists  important  geometric  and  loading  parameters  for  the 
75  surface  crack  cases  and  Table  A-2  lists  similar  information  for  the 
85  corner  crack  cases.  Table  A-3  lists  10  coefficients  for  each  indi¬ 
vidual  case.  These  coefficients  result  from  the  original  finite  element- 
alternating  method  calculations  and  must  be  used  as  input  data  for  the 
computer  programs  described  in  Appendices  B  and  C.  A  list  of  coefficients 
for  cases  different  than  those  considered  in  the  present  report  is  con¬ 
tained  in  Reference  10.  These  coefficients  can  also  be  used  in  the 
computer  programs  of  the  following  Appendices. 

The  actual  plate  modeled  by  the  finite  element  -  alternating  method 
had  a  thickness  of  1.0  and  a  radius  0.5.  The  objective  was  to  compute 

stress  intensity  factors  for  crack  pressures  p(x/R)  *  A  (x/R)n  = 

Ar  n 

— -  x  .  Since  R  =  .5,  the  actual  loadings  considered  for  the  current 

R 

model  were  p(x)  =  (An/Rn)xn  »  1,  2x,  4x2,  8x3  for  n  »  0,  1,  2,  3.  In  a 

few  cases  denoted  in  Table  A-2,  the  pressure  was  p(x)  =  1,  x,  x2,  x3. 

These  results  were  converted  to  the  standard  loading  by  multiplication 

by  an  appropriate  constant.  Thus,  all  of  the  results  given  in  Tables 

4  and  7  are  for  p(x/R)  =  1,  (x/R) ,  (x/R)2,  (x/R)3  which  is  equivalent 
2  3 

to  p(x)  =  1,  2x,  4x  ,  8x  for  the  actual  model  (R  =  0.5). 
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Table  A-l.  Surface  Crack  Cases 


Table  A-l.  (cont'd) 


Case  Crack 

Number  Type  a/c  2a/ 1 

Loading  Case,  p(x)  = 

2  3 

Remote  Tension  1  2x  4x  8x 

150  Surface  1.5  .2 

X 

151 

X 

152 

X 

153 

X 

154 

X 

155  .4 

X 

156 

X 

157 

X 

158 

X 

159 

X 

160  .6 

X 

161 

X 

162 

X 

163 

X 

164 

X 

165  .8 

X 

166 

X 

167 

X 

168 

X 

169 

X 

17C  1.0 

X 

171 

X 

172 

X 

173 

X 

174  x 
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Table  A-2.  Corner  Crack  Cases 


Case  Crack 

Number  Type  a/c  a/t 

Loading  Case,  p(x)  = 

2  3 

Remote  Tension  1  2x  4x  8x 

175  Corner  .5  .2 

X 

176 

X 

177 

X 

178 

X 

179 

X 

180  .4 

X 

181 

X 

182 

X 

183 

X 

184 

X 

185  .5 

X 

186 

X 

187 

x* 

188 

X* 

189 

X* 

190  .9  .2 

X 

191 

X 

192 

X 

193 

X 

194 

X 

195  .3 

X 

196  .4 

X 

197  .5 

X 

198  .6 

X 

199  .7 

X 

2  3 

*  p(x)  for  these  cases  was  x,  x  and  x 

52 

.  4 


Loading  Case, 


Remote  Tension  1  2x 


4x2  8x3 


Corner 


*  p(x)  for  these  cases  was  x,  x  and  x 


/ 


Case  Crack 

Number  Type 


a/c 


a/t 


Loading  Case,  p(x)  s 


Remote  Tension 


2x 


4x 


8x" 


224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 


Corner  .  7 

.8 


.9 


1.5  .2 


.4 


.6 


.8 


x 


x 


X 


X 


X 


x 


x 


x 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 
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Table  A-3.  (cont’d) 


Table  A-3.  (cont'd) 


Table  A-3.  (cont'd) 


r 


i 


! 

b. 


60 


Table  A-3.  (cont'd) 


Table  A-3.  (cont'd) 


APPENDIX  B 


MODE-ONE  STRESS  INTENSITY  FACTOR  CALCULATIONS 

1.  INTRODUCTION 

Tables  4  and  5  and  Tables  7  and  8  In  an  earlier  section  of  this  report 
list  dimensionless  stress  intensity  magnification  factors  for  corner  and 
surface  cracks,  respectively,  at  holes  in  plates.  These  magnification 
factors  are  shown  in  the  interest  of  brevity,  for  30° increments  along 
the  crack  border.  It  is  entirely  possible  that  at  some  future  date  simi¬ 
lar  information  for  other  locations  might  be  needed.  This  appendix  de¬ 
scribes  two  computer  programs,  one  written  in  FORTRAN  and  the  other  in 
BASIC,  which  generate  mode-one  stress  intensity  factors,  Kj,  at  one 
degree  increments  of  parametric  angle  along  the  crack  border. 

2.  PROGRAM  DESCRIPTION 

Programs  KFRTN  and  KBSC,  whose  listings  are  found  in  Tables  B-l  and 
B-2  respectively,  perform  identical  functions.  The  difference,  of  course, 
is  the  computer  language  used.  Program  KFRTN  is  in  FORTRAN  and  KBSC  is 
in  BASIC  and  this  latter  program  was  written  and  tested  on  the  Hewlett- 
Packard  series  of  desktop  computers.  Both  programs  are  structured  in 
a  standard  fashion.  The  first  section  reads  and  writes  input  parameters. 
The  second  section  determines  the  parametric  angles  for  which  K^.  will 
be  calculated  along  the  crack  border  depending  upon  the  crack  location 
and  shape.  The  third  and  final  section  performs  the  lengthy  algebraic 
calculation  to  arrive  at  K^. 

3.  PROGRAM  INPUT 

Table  B-3  lists  required  variable  inputs  determined  for  the  particular 
case  being  considered.  Figure  B-l  graphically  depicts  these  parameters 
in  more  detail.  The  authors  suggest  that  any  user  of  these  programs 
exercise  special  care  to  be  absolutely  certain  these  parameters  are 
correctly  input  (sketches  similar  to  those  of  Figure  B-l  are  particularly 
helpful) . 
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The  user  should  also  be  aware  of  certain  values  fixed  in  the  software 


which  must  remain  as  shown.  All  calculations  are  for  a  Poissons  Ratio  of 
0.25  and  a  Modulus  of  Rigidity  of  12x10^  psi.  Stress  intensity  factors 
obtained  from  the  stress  function  coefficients  are  given  in  Table  A-3 
for  specific  loadings.  For  the  remote  tension  cases,  the  applied  load 
was  1  psi.  For  the  crack  face  pressure  problems,  the  applied  load  is 
p(x/R)  =  Aq,  A^(x/R) ,  A2 (x/R)^,  Aj(x/R)4,  where  AQ  =  A.^  *  A2  *  A^  ■  1 
and  R  =  0.5.  Other  loadings  and  hole  diameters  can  be  considered  by 
scaling  the  resulting  stress  intensity  factors  in  an  appropriate  manner 
(multiply  by  desired  remote  stress,  pressure  load  Ar,  and/or  relate  crack 
dimensions  through  a/c  and  a/t).  In  all  cases,  one  must  maintain  2R  «•  T. 

4.  PROGRAM  OUTPUT 

Table  B-4  shows  a  sample  output  of  Program  KBSC  (output  of  KFRTN  is 

similar).  The  mode-one  stress  intensity  factor,  K^,  is  printed  for 

increments  of  parametric  angle,  where  this  angle  is  always  measured 

from  the  crack  semi-minor  axis  as  shown  in  Figures  5  and  6.  For  the 

o 

corner  crack,  will  be  calculated  for  0  to  90  in  one  degree  incre¬ 
ments  but  for  the  surface  crack  will  be  calculated  for  -90°  to  +90° 

(all  along  the  crack  border)  in  two  degree  increments.  Note  that  in  the 
outputs  of  programs  KBSC  and  KFRTN,  the  parametric  angle  $  is  always 
measured  from  the  semi-minor  axis,  while  in  Tables  4  and  5,  $  is  always  de¬ 
fined  as  the  distance  from  the  front  surface  for  the  corner  cracks  as  in 
Fig.  5.  Thus,  for  a/c  >  1.0,  both  the  tabulated  and  computer  program  output 
have  the  same  definition  for  while  for  a/c  <  1.0,  the  angles  are  com¬ 
plementary.  Since  a/c  always  exceeds  1.0  for  the  semielliptical  cracks 
considered,  there  is  no  inconsistency  in  the  two  definitions  for 


65 


Table  B-l  Program  KFRTN  Listing 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


THIS  PROGRAM  CALCULATES  MODE  I  STRESS  INTENSITY  FACTORS  AT  ONE 

DEGREE  INCREMENTS  Of  ELLIPTICAL  ANGLE  ALONG  THE  BORDER  Of  A 

SEN  1  ELL 1PT1 C AL  CRACK  INTERSECTING  A  FASTENER  HOLE  IN  A  FLAT, PLATE 


INPUTS  ARC 

NAXIS  *  CRACK  ORIENTATIONCSEE  EXPLANATORY  FIGURES) 

A  =  CRACK  SEMI-MAJOR  AXIS 
B  =  CRACK  SEMI-MINOR  AXIS 

D1  =  DISTANCE  FROM  PLATE  CENTER  TO  CRACK  CENTERCSEE  FIGURES) 
CCD  =  STRESS  FUNCTION  COEFFICIENTS*  INPUT  IN  ORDER  OF  TABLE 


OUTPUTS  ARE 

ANGLE  =  ELLIPTICAL  OR  PARAMETRIC  ANGLE  FROM  CRACK  MINOR  AXIS 
KJ  =  MODE  ONE  STRESS  1NTENSITT  FACTOR  FOR  ONE  PS1  REMOTE  STRESS 


DIMENSION  CMO) 

G=  12000D0D. 

PI  =  3  .1*159 

T1  =  1  .r* 

RE  A  P  (5  /  '  )  NAXIS/A*B*D1 
WRITFC6#?)  NAXIS 
WRlTr(6*3)  A 
WRI1EC6/*)  B 
VRITE<6/5>  D1 
RE  AO  (5  #6)  tC(II*l=1,8> 

READC5.6)  (C(I)*I=9*10) 

WRITE <6*I1) 

NR1TEC6.7)  <C(I)*1=1*10> 

WRITE (6,8) 

DETERMINE  ANGLES  FOR  WHICH  STRESS  INTENSITY  FACTORS  WILL  BE 
CALCULATED. 


0PP=.5«TT«D1 
OP P= AM A XT COPP/-A) 

CP  P=  AM  INI  (OPP/A) 

T1  =  ASINCOPP  FA) 

OPP=.5*TT-Dl 
OrrrAMAXl CCPP/ -A ) 

OPP= AMIN-  (OPP/A) 

T?=ASINCOPP/A) 

DT=(T2«T1)*2./Pl 
71=11*180. /PI 
If  CHAXIS-2U0/20/30 
3C  PRINT  1  •"00/  NAXIS 
1C~C  FORMAT (IX/'NAXIS  NOT  T  OR  2**J10  ) 
2  C  CONTINUE 
T 1 =T  T -90. 


GO  TO  AO 
IT  CONTINUE 
*r  CONTINUE 


DO  ir.n  1=1/91 
ANG=(I-1)*DT-T1 
CAL  C  =  -AIJG«?9. 

1N=CALC«PI/U'. 

AKI  =  f .*GFf A*P)«S0RT(P1/CA*B>)»  CA»A* (SIN  CTH)  )*«?* 

1B«B« (COSCTH) >♦•?)««. 25* (CCD 4C (3 )«SIN(TH)/B-*.« 
2C(A)*(COS(TH>)**2/<A«A)-A.«C(6>*(SIN(TK) )«•?/ (B«B) 

3-A.*C(B)*(C05«TH))*»?*SIN(TH)/CA*A«B)-*.*CllO)* 

IN  (TH))**3/8«*3AC(2 >*C05(7H>/A  +  C<5)*.C0S<TH)*SIN(TH>/CA*B}~4.* 
SCC7)*(C0S  (TH>/A)»*3-*.«CC9>*  CCOSCTK) /A)*CS1N(TH>/B)*«2> 

"R 1TF (6/9/  ANG/AK 1 
ANG=r..n 
CALC -0.0 
AKl=r.O 
TH=?.P 
17T  CONTINUE 

1  FORMAT  (1  1C/3M0.*) 

5  FORMATCIHC/’NAXIS  *• /12/F) 

FORMATCIHC# 'CRACK  SEMI-MAJOR  AXIS  s'*E5.2*/> 

FORMAT  ClHC/’ CRACK  SEMI-MINOR  AXIS  *,/E5.2*/) 
fOR«AT(lHr/*DISTANCE  TO  CRACK  CENTER*  01  *'*F$.2*//) 

FORMAT  (I0EI3.O 

FORMAT  C//*6X*'  ANGLE  '*6X*'Kl»*/> 

.  FORMAT  crT0.T*FTr-.'3) 

11  ^gpATOMC',sms5  FUNCTION  COEFFICIENTS*  CC1)  TO  C(10)'*/> 

END 


3 

* 

5 

6 
7 
t 
9 


V 


Table  B-2.  Program  KBSC  Listing 


1 0 

DIM  Hi 

5  3 .  BI 

5  J-II 

4  3.  F  t  4  3  <  Cl  10  3.  IK  10  1 

t‘l  1 

!  'P  l  NT 

"  T  H 1 S 

PPOG 

RRM  CRLCULhTES  MODE 

ONE  STRESS 

IN TENS 

>  I T 

:ii±t 

1  OR  1  - 

1  TO 

10 

40 

or  1 3=0 

00 

mi:  :r  i 

c  o 

MOP  •• 

]  NPM 

ORRC 

K  SEMI- MR JOE  RND  ' -1 

Ml  MINOR  AX 

ES : H . B ' 

n. 

70 

INPUT 

R .  B 

Bit 

pr.  in  r 

L1H.'. 

"  SEN  I 

. MAJOR  AXIS  =">0? 

■  El'll -MINOR  R 

FIS  ="! 

!  B 

00 

HI  1  1  0 

I  m  s  i  II  2  3:-;  0 .  1 

I I  ii  wr  3  i--0. 0] 
i  :•!!  kc  4  ]-0. 2 
120  HI  5  3-0.  1  1 
1  40  ML  6  3-6. 02: 

130  I  lf  7  ]=0. 3 

5  00  HI  0  3-0. 21 
1  V ill  III  9  3-0 . 12 
J  00  Hi:  10  I--0.03 

1  •  'i:i  lij "  1  HP  U[  0.1  RE  S3  FUHCTIOH  COEFFICIENTS" 

2uu  Hr  ill  3000 

flu  inn')  L  1112'.  "|  00  TO  003 "  .  L  1  H  I 
' . ' u  IMF  1  =  1  TO  10 

2  ;n  i  IF. I'll  2 

:  in  iu  of  ■'(  ' .  Hi  1  i. 

.  "-U  I  HIJ  DR FI.i 

2r  O  I  III 'in  cu:i 

70  i i  HP]  7 £'  (2>  2 SO) i  l  1  J 

2::u  i  ij FLIRT  El  3. 3 

2”  4 1  I  IEFT  1 

4O0  iM'.P  "EOF'  CORRECTIONS  ENTER  1.  II  HO  CORRECTIONS  ENTER  0"5 

MO  INF IJT  K 

420  IF  I  =1  THEN  340 


f  LIN 


04 1  got  o  Too 

040  lm.p  "ENTER  O0EFI  1CIEHT  NUMBER,  i  1 FF )  =? .  AND  PRESS  CONT-EXECUTE 
SCO  l  Of 

01.0  l.i  ISP  "NR  MS  " ; 

370  INPUT  H2 

440  i;  12000000 

:•  «0t  Hi  OP  "Hi STANCE  TO  CROCK  CENTER  =  "? 

400  INPUT  M 


410  1-1 
4  20  PE: IN  I 
4  70  PR  INI 
440  PRIM) 
4  00  PE' 1IIT 
4i.ii  I  P  INI 
4  70  PE' I  HI 
•  i..0  Pi-  1111 


L  illi  .  "IlHFlb  =  "  •  M2 
"  G  -  "  5  0 

"  I'  I  SI  ONCE  TO  THE  LETU  P  CENTER  "  5  D 1 

"PL ME  TH)  Cl  HESS  =  "  5  T 

LIII4.  "tlNGI  I  =  P  HR  RITE  TR  I  P  llNi.,1.  \.  I  ROM  MINOR  RFISMLIHl 
”P1  MODE  -ONE  S  TRE':-S  II 1 1 II P.  i  I  Y  FACTOR "  i  L  INI 
LIN.  .  "llUGL  E "  j  SF’OSJ  "  I  J " 1  IN.M 
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Table  B-2.  (cont'd) 


4 

90 

0 1  = 

0 . 5 

*t+m 

5 

0Ti 

110  = 

-ft 

5 

10 

IF 

0 1  < 

ft  9 

THEN 

530 

cr 

.j 

20 

GOTO  5 

40 

cr 

30 

0 1  = 

ft  9 

cr 

•J 

4  0 

IF 

o  l  > 

ft 

THEM 

560 

c 

50 

iicn 

ft  5 

70 

S  b  n 

01  = 

II 

5 

7  U 

f  1 2  = 

0 1  / 

ft 

c: 

30 

IF 

02  = 

1 

THEN 

600 

l: 

90 

Gu  1 

0  6 

20 

b 

00 

i  1  = 

pi 

2 

6 

1  0 

till  1 

0  r 

i  3  0 

b 

2  u 

T  1  = 

■ft  i  N 

<0 

2  •  i,iP 

'll . 0212 

f. 

0 

1 1 1  = 

■  D .  51 

'*  T 

. ]i  1 

b 

40 

1  F 

0  1  < 

ft  9 

THEN 

660 

b 

60 

’,0  1 

ij  * 

.  7  0 

b  Ci 

1 1  ;i  -■ 

ft  9 

b 

70 

i  i 

o  i : 

ft 

1  HEN 

690 

b 

30 

GUT 

0  7 

00 

b 

•  M 

I  1  1  :• 

ft 

1  lit 

1  .! 

0  1  ■ 

'ft 

1 1  0 

11 

l.,l  J .. 

1 

THEN 

.“30 

, ■ ■  m  i, ii  10  750 
,"ju  i ::  pi.  2 


'4  0 

601  0  760 

750 

12  ft  T  N  <  Q2-  S 

OP 

<1-02T2 

.1  1 

4,0 

l.i 3~ (  T ''+ti  , 

2  ■ 

I 

770 

11=11*180  P 

1 

,3.11.1 

IF  112  =  2  THE 

N 

300 

;"‘9  i'.i 

6OT0  310 

iitn.i 

ri=Tl-90 

1  0 

TOP  1  =  1  Tli 

91 

320 

ft  3= < I  -  1 > *D2 

...  T 

1 

3  ;0 

ti 4  —  -  R3  +  90 

0411 

ft  5:  II 4*  P  I  1  8 

0 

rib  IJ 

l"  3-8  *  G  <  ft  B 

i  ft 

ftp '  p  I . 

■  M  '  [ 

j  .'  *  <  ft  •  1 1 

*• '..  ft  III' 

ft  5  > 

1 2 + B  *  B  *  <  C  0  S  C  ft  5  T 

)  12 )  10 

:::bn 

1  4  =  Cl'  1  ]  +  !..  [  :• 

i  ■ 

1  I  I  '  ft  5 :. 

.  B . 

THI  4  J<  ( 

30 SC  ft! 

i  '  >  t  c 

:VCft*fi) 

:V,"i  i 

r  5-~4*C[  6  1* 

(  '• , 

i  h  <  ft  5  ■ :. 

•r  2 . 

B*ft;.i~4  * 

218  J* 

.COS' 

ft  5  )  >  1 2  *  S I H  <  ft  5  >  7 

c  ft*  ft*  c 

.  .11 

i  -7-'-  4  *  cr  io  i 

t  t 

3 1 N  i  ft  5 , 

■  13. 

n  1 3+cr  = 

3*  COS 

1 15) 

'ft 

"i 

i  7-  CL  5  3*1  1!': 

,  h 

■*  SIN ■ 

H5>. 

'.  f1*B  1 . 4 

*|'.[  7  ] 

-  ■  ft  Oi 

5  C 115)  •■■'ft)  13 

9  MU 

1  4 ---4  *C[ '4  J  ► 

(  I. 

OS  <  h5). 

ft  >  * 

1  H  •'  ft  5  ' 

•B  '  12 

9  1  i  1 

i  •  ■  1  3*  ■:  14  1 1 

J  ‘ 

1  6  +  F  7  + 1 

;  j 

97U 

HP  HE  '"2-  'l 

0 

II 3  <  K9 

'  •  M 

i  iiPllftT  F 4 .  O 

>  F 

1  2 .  4 

740 

Ml  7.1  1 

.0 

1  NH 
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Table  B-3  Input  List 


Input 

Description 

Determined  From 

NAXIS 

Define  the  crack  orientation 
vlth  respect  to  the  plate 

Figure  B-l 

A 

Crack  semi-major  axis 

Table  A-l  or  A-2 

B 

Crack  semi-minor  axis 

Table  A-l  or  A-2 

D1 

Locates  the  crack  center 
with  respect  to  the  plate 

Figure  B-l 

C(I) 

Coefficients  to  calculate  Kj  where 

C(l)  -  COO,  C(10)  -  C03 

Table  A- 3 
(reading  from 
left  to  right) 

Corner  Crack  (a/c  £  0.9,  2R  =  1 


Table  B-4.  Sample  Output,  Program  KBSC 


nr.  rruGKiii-i  uiu  ulhtes  rii iji!  n m:  .  1 1’ ■  .  .  mu: ns i iv  f fit t of 


mi  in i .top  h: : i s  i  SLi'i i - 1 1 1 1 10 r  o : : ; •  o . 5 


i in  Hi  CO 3 

/0 ,  0 . fv.i 

•  ssct-i  u 
.  1 37E-  1  0 
■I.  :4  RE  - 1 1 1 
■  ■of:-'  1  II 

.  V;  1 1  :  -  l  I 
FOEi- 1 1 1 
i . lObE-10 

25?E~  1  S 

' .  1  0  ,■  E  -  1  1 

i SHUOOOU 

: • . 1 1 1  si  t:  to  in  •.  knot  ce nil r  m' 

i  mi  mil  KNiv  i 


■  I  MR  Ml  ii  I  !■'  I  .  IIHnLI.  f  I-- 1. '  I'i  11  1. 1 1 IJ  I’  I  i,  .  i 
illjli,  ■  OUR  '  .  I  l-'l  I  Hit:  I  Ml  I  i  1  HI  10!' 


APPENDIX  C 


CRACK  OPENING  DISPLACEMENT  CALCULATIONS 


1.  INTRODUCTION 

Crack  opening  displacements  can  be  easily  determined  using  input 
similar  to  those  for  the  mode-one  stress  intensity  factor  calculations 
of  Appendix  B.  Reference  10  fully  describes  the  formulation  of  a 
lengthy  yet  simple  algebraic  expression  which  results  in  crack  opening 
displacements.  This  appendix  describes  a  computer  program  which  will 
generate  such  displacements  at  certain  locations  on  the  crack  plane. 

2.  PROGRAM  DESCRIPTION  AND  INPUT 

Computer  program  COD  (listed  in  Table  C-l)  is  written  in  BASIC  on 
a  Hewlett-Packard  desktop  computer.  The  derivation  of  the  algebraic 
expression  which  yields  crack  opening  displacement  includes  A,  B,  and  C(I) 
which  are  exactly  as  described  in  Appendix  B  for  Programs  KFRTN  and 
KBSC. 

3.  PROGRAM  OUTPUT 

Program  COD  output  includes  crack  opening  displacements  calculated 
and  printed  for  specific  points  along  rays  at  10°  increments  of  para¬ 
metric  angle,  $,  over  90°  of  the  crack  where  $  is  measured  from  the 
crack  major  axis.  Figure  C-l  depicts  a  few  of  these  points  for  clarity. 
Each  ray,  whose  length  might  be  called  r^,  has  output  points 

O.Or^,  O.lr^,  ....  0.9r^,  0.95r^,  0.98^  . 

These  points  coincide  with  the  description  required  for  Benchmark 
Problems  and  an  example  can  be  found  in  Reference  3. 

Table  C-2  shows  a  sample  output  from  Program  COD.  The  coordinates 
of  each  point  are  described  in  terms  of  x  and  y  (Figure  C-l)  rather 
than  r^.  As  with  the  Appendix  B  calculations,  for  the  remote  loading 
cases  the  uniform  remote  tension  was  one  psi. 
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Table  C-l.  Program  COD  Listing 


1  O  D I M  XI  12],  YI  1  2  ] ,  N  [12],  Cl  1  0  1 ,  El  1 0  J ,  111  1  0  ] 

2u  PRINT  "THIS  PROGRAM  CALCULATES  CRnCR  OPENING  DISPLACEMENTS" j LIN2 
30  FOR  1=1  TO  10 
40  El  I  1-0 
5U  NEXT  I 

Co  D I SP  "  I  HP  UT  C  ■  ACK  SEM I  -MAJOR  AND  SEN  I  -  M I  NOR  AYES :  A ,  B "  5 
70  INPUT  A, B 

CO  PRINT  LINE, "SEM I -  MAJOR  AXIS  =" 5  A? "SEM I -MI NOR  AXIS  ="!B 
00  U-0.23 

100  PRINT  LI  NX, "POISSON’S  RATIO  =";U5LIN2 

110  NC 1 3=0 

120  HI  2  1=0.1 

ISO  HI.  3  3=0. 01 

140  UC  4  3=0. 2 

150  NI 5  ]=0. 1 1 

ICO  111  6  1=0.02 

170  HI  7  1  =  0.3 

ISO  NIS  3=0. 21 

1.90  NC  9  3=0.12 

200  NI  101=0.03 

210  H ISP  "INPUT  STRESS  FUNCTION  COEFFICIENTS" 

220  WAIT  3000 

2 30  PRINT  LIN2,"CO0  TO  COS",  1  INI 
240  FOR  1=1  TO  10 
250  FIXED  2 

2C0  It  ISP  "  C  "  ,  NI  I  1,  "  =  "  5 
270  STANDARD 
230  INPUT  Ct I J 
290  WRITE  (2, 3U0  ' CL  II 
3O0  FORMAT  El 3. 4 
31.0  NEXT  I 

320  I.USP  "FOR  CORRECTIONS  ENTER  1,  II  NO  CORRECTIONS  ENTER  0"5 

330  INPUT  K 

340  IF  K = 1  THEN  SCO 

3 50  GOTO  380 

SCO  D J  SP  "ENTER  COEFFICIENT  NUMBER,  f.  .  XX>=?,  AND  PRESS  CONT-EXECUTE" 
370  STOP 

300  01  =  <  A 1 2 ) *  1 B  t  2  > 

390  A 2=  <  A 1 2  >  *  I  fit 4  ) 

400  A3- ( Ht4>*>  BtC> 

4  1  O  A4  =  (  A  't  2> +  •:  B  t  C> 

4  20  05=  (.  A 1 4  >  •?  B't'8 ) 

4  ;0  h 6--'.:  A 1 2  >  =  BtS  "> 

440  07  =  (At  4  ■  +  ■  Bt2> 

450  A 3= < A  t  4 ) * 1 B  t  4  :> 

4 1 .0  N9-  <  A t  6  ">  * '  .  BTC  > 
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Table  C~l.  (cont'd) 


•t  :V'  i  ;  I  •  i'i  '■  r.  .•  •  •  i  I  -i  ' 

i:  i  : . '  1  f  |  .  1 

•;  'it  r  •  h ' C-  •••)•. U'  • 

' -I  Ml  :  ;  ■  ti  rf;  .  ■  .  u  T  I 

!  1 1  !  Ilf.'  1  \  '  in 

I  : '  1  i  n  -  •.  1  1  •>»  1  I'  ll 

i- :tn  1. 1 !  4.  i  “i'i  !!•'!  i!  n  i!  f  iii!1,!  i  i  !•  1 1  h  :n,  i.  •).-  M.i:,  -  .  j .  j  'Me; 

‘.•i.:  ,  !!•'  .11  1  1 

!  I  if.  1 

'  1  i  i  I.’  1  ,  14  i -tii  i 

n  !  :  1-11  I  111  M  njll 

;  I  i  j  !  !  il!  II  I 


•Ml  j 

Mil  m.i 

un  j 

i  ■  U  .  M05 

1U  i.i 

ii|  0  t.  30 

•  k  i  ; 

L  lit .  M 

,  .Ti  i 

1-  .Wh*M 

mM  i 

'=■  >[:<!'  1 

Ml 

1  i  VC  Hu 

♦  l» 

S  i  m;  hr i 

, ’l.l  i., 

■  1  1  '.t  ,'lfl 

:  *  1 1  r, 

nlii  V'U-i 

•  ■ :  i  . ' 

i  .  n 

'Hi' 

,  i  .l-ii 

:  *  i  * 

i  •  hi:'.  ■  .  n 

mi 

.  'll': 

■ '  II  i  I  1  ]  *  •,  •)  Ill  I  !  I  II.'  Ill  2,  ,i  Hit  6  J 

.  .,ii  I  i  l  10  )  •  !  >71  l.l'',/  lY.-H.l  \w}-  I."','!  .1  11/ "'il  ,"f|i:-.-C.r  2IH4*HC  J] 

.  ‘ .  n  '■  1 1. 1 .  .1*/ > '  1  H t  l  -i  I  *  11''  '  in  t  I  -M  ■ : :  •  ■ . ’  [>  1  *  Cf  4  3*  1 6*  <  X  [ 

I-  It  'j*;:  ,1  i  1 1 '  1 1 1.  [  :  J  *  ]  i  ■  ■  1  ,  I  Ml,  ..VI  -B1 . C[  7  1 *24  ■>:  >! 

,  i  I  .  I  .  1*1’  ■  . :  I.  I  !■', 

'  ;  •  >’  i  .  •  U  ,  '  '  I'H  !  '  I  0  I  M  V  ' 

■  n  !i:  i  i 

:  n,  1'fUU  L  Ui  I  !  Ii  l  t  <  ' !/  "  5  f'H  In!  "V  -  lit.:;  1 1  I  ‘ . I  ' I  "  j  I  .  Hll 

/  :  n  ui.'  T  l  l  ■  -  m  u  ; ;t  i  /  I.  i  m  /  ] -  Ilf  1  /  . 

:  .  >  ■  II:  :  /  1-  u 
: '  .  i  i  ,iK  l-l  in  11 


Table  C-2.  Sample  Output,  Program  COD 


I  HIS  I'i-'t iijf'Hl'1  CHLCULh  I  ES  CkFlCI  OPEN  I  HR  IlISPI 

SEMI  HU.  I  OR  RXIS  -  1  SEMI-  HIM  Ok  P  ,!'•  =  R 

Pu  i  •  itr  PR  I  ]  II  --  R .  .  '• 


MU  1  M 

1. 1,1  J 

•  ‘ 

T7RE 

-89 

6  -  2. 

■ji.'.ME 

- 1 R 

2.  1 

SERE 

1  R 

0 

47RE 

-  1 0 

- 

90RE 

-  1  R 

;> » 

•:  1  RE 

-11 

-  2  .  2 

.'Eire 

-lit 

1  -  1 

i.  RE 

-  1  R 

«  2 

5  o' RE 

-  1  1 

1. 

R  7'u[' 

-  1  1 

I'liPtiHi  I  !•' 1 1  ItHl.EE  EkOI'1  HfURk  I  M  !  1  n  Dt 


D  l  SI  E 


y  „  00  0 

R .  ti0i  i 

4.oi  64  e: 

11 

it .  1  RR 

R.  RRO 

4.761  IE 

0 

0 . 20M 

u ,  RRR 

4 . 6 1 63 E 

R 

R.  JRR 

y  I,  0  M  u 

.  soyoL 

R 

R.  4  Mil 

0 . 00  M 

•i  .  R:-:6RE 

0 

y  i .  0  U  0 

R .  0Ro 

4.71  4  EE 

0 

R .  60R 

0 . 000 

: . 27 75E 

Vi 

0. 7R R 

R .  RRR 

2. 7  7.., 9 E 

0 

0 . 3RR 

0 . 001.1 

. 'mrCE 

0 

Li .  'TOR 

0 . 000 

1  .  -5 ' <  r 

0 

h.  950 

11 . 001  1 

!  .  R  :::E,F 

11 

i  i ,  9SR 

0  „  000 

4  0.  "SIS 

0 

75 


_ ft CEMENT 


I 


Table  C-2.  (cont'd) 


PARAMETRIC  ANGLE  FROM  MAJOR  AX I: 


10 


DEGREES 


DISPL 


000 

0 . 000 

4.8164E-08 

098 

0 . 009 

4 . 7697  E  -■  0  8 

197 

0.017 

4 . 63 36 E -08 

£95 

0. 026 

4. 4134 E -08 

394 

0.035 

4.  1 150E--08 

492 

0. 043 

3. 74 36 E- 08 

591 

0.  052 

3 ,  3  0  8  5  E  “  0  8 

689 

0.  061 

2 . 7  9  4  9  E  0  8 

7* !-!  I~l 

0 . 069 

2 . 2075E-08 

886 

0.078 

1 . 495 IE -'08 

9  3 6 

0. 082 

1 . 0308E-08 

9  6  5 

0, 085 

6 . 4 1 :  3  4  E  -  0  9 

• 

• 

• 

• 

• 

• 

m 

• 

• 

■  ANGLE 

FROM  MAJOR 

AXIS  =  90 

DISF'L 


0 . 000 

0. 000 

4.81 64E--0 

0.  000 

0. 050 

4.8380E-0 

0 . 000 

0  .  1  00 

4 . 8082E-9 

0 . 000 

0.  150 

4.7235E-0 

0. 000 

0.200 

4 . 5783E-0 

0 . 000 

0.250 

4,36  3 9  E — 0 

0 . 000 

0. 30O 

4 . 0  6  6  2  E  —  O 

0. 000 

0. 350 

3 . 66 1 4E--0 

Cl.  000 

0. 400 

3. 1Q3 IE- 6 

0. 000 

0.450 

2.2743E-8 

0. 000 

0.475 

1.6365E-0 

0 . 000 

0 . 490 

1 . 945SE-0 
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Figure  C- 


.  Program  COD  Output  Points 
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